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ABSTRACT 

Nonequilibrium (timc-dcpcndcnt) cooling rates and ionization state calculations are 
\^ ' presented for low-density gas enriched with heavy elements (metals) and photoionized 

\ by external ultraviolet /X-ray radiation. We consider a wide range of gas densities and 

metallicities and also two types of external radiation field: a power-law and the extra- 
galactic background spectra. We have found that both cooling efficiencies and ionic 
composition of enriched photoionized gas depend significantly on the gas metallicity 
O ' and density, the flux amplitude, and the shape of ionizing radiation spectrum. The 

cooling rates and ionic composition of gas in nonequilibrium photoionization models 
. differ strongly (by a factor of several) from those in photoequilibrium due to overi- 

onization of the ionic states in the nonequilibrium case. The difference is maximal 
at low values of the ionization parameter and similar in magnitude to that between 
£NJ ■ the equlibrium and nonequilibrium cooling rates in the collisionally controlled gas. 

J> ' In general, the nonequilibrium effects arc notable at T ;$ I0 6 K. In this temperature 

, range, the mismatch of the ionic states and their ratios between the photoequilibrium 

C^l ■ and the photo-nonequilibrium models reach a factor of several. The net result is that 

the time-dependent energy losses due to each chemical element (i.e. the contributions 
to the total cooling rate) differ singificantly from the photoequilibrium ones. We ad- 
0^ ' vocate the use of nonequilibrium cooling rates and ionic states for gas with near-solar 

(and above) metallicity exposed to an arbitrary ionizing radiation flux. We provide a 
parameter space (in terms of temperature, density, metallicity and ionizing radiation 
flux), where the nonequilibrium cooling rates are to be used. More quantitatively, the 
nonequilibrium collisional cooling rates and ionization states are a better choice for 
the ionization parameter log U ^ — 5. The difference between the photoequlibrium 
' and the photo-nonequilibrium decreases with the ionization parameter growth, and 

the photoequilibrium can be used for ionization parameter as high as log U Si — 2 
for Z ^ 10~ 2 Zq and log U ^ for Z ~ Zq. Thus, the nonequilibrium calculations 
should be used for the ionization parameter range between the above-mentioned val- 
ues. In general, where the physical conditions are favour to collisional ionization, the 
nonequilibrium (photo) ionization calculations should be conducted. 
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1 INTRODUCTION 



Recent observations of heavy elements (metals) in 
the intergalactic medium (IGM) give rich information 
about star formation and enrichment history of the Uni- 



verse |Madau et al. 


1 19961; Pettinil 19991: Schave et al. 1 20031: 


Nicastro et al. 1 2006 


: Simcoe et al. 200a Danforth & Shulll 


20081). In order to understand these observations, it is nec- 
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essary to know which processes are responsible for the ob- 
served ionic composition and thermal state of the IGM. 
In the absence of ionizing background radiation, these 
processes are controlled by collisions and determined by 
gas metallicity and temperature. However, when galax- 
ies an d quasars form they p r oduce strong ionizing pho- 
tons (lHaardt fc Madaul Il99fj , l200lt iMiniati et al~l 12004 
iFaucher-Giguere et al. I l2009l b Therefore photoionization 
and photoheating are to be taken into account. 

Calculations of the cooling rates of astrophysical plasma 
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in the collisional ionization equilibr i um ( C IE) were 



form ed by many autho rs dHousel 1964 1_ Cox fc Tucke 



1969 



1982 



Raymon d et al. I 1976 
Gaetz fc Salpeterl [l98S 



Shull fc van Steenbere 



Bohringer fc Henslerl 



19891; ISutherland fc Dopital 1 19931; lLandi fc Landinil 1 1991 



Benjamin et al. I l200ll ; iBrvans et al. I \2QQ$ \. However, 



calculations of the time-dependent ionization of metals and 



from 


the CIE states (iKafatosI Il973; IShaDiro & Moore 


197£ 


; Edgar & Chevalier! 1986; 


Schmutzlcr & Tscharnutcr 


1993 


; ISutherland & Dooita 1993 


; Gnat & Sternberg 2007h. 



In some conditions, ionizing radiation can be important 
for the ther mal and ion i zation states of primordial and 
enriched gas. lEfstathiou I (|l992h showed that the presence 
of UV radiation can strongly suppress the cooling rate of 
primordial gas. In th e photoionization equilibrium case, 
W lcrsma et al. I (|2009l ) came to a similar conclusion for 
enriched gas exposed to the external ionizing radiation. It is 
obvious that ionizing radiation forces the ionic composition 
of gas to settle on to photoequilibrium. However, during 
transition from a pure collisional to a photoequilibrium 
case the nonequilibrium effects are expected to play a role. 

In general, numerical simulations of metal ionization 
states in the IGM should take into account the effects of 
ionizing background. However, this requires tracing huge 
number of metal species and their ionization states, as well 
as complex calculations of the cooling and heating rates of 
the enriched photoionized gas. Therefore, simulations often 
consider either a limited number of metal ionization states 
or cooling rates o f collisionally ionized gas. For instance, 
iRaga et al. 1 (|l997f ) studied a radiativ e bow shock using a 
limite d ionization and cooling network; lYoshikawa fc Sasakil 
(|2006h used the CIE cooling rates to study the ionic compo- 
sition of a warm-hot intergalactic medium. iHughes fc Singh] 
|l994h presented a comprehensive nonequilibrium ionization 
analysis of the supernova remnant. Recently, self-consistent 
calculations were success fully performed in one-dimension. 
iGnat fc Sternberg] (|2004T l studied the metals photoioniza- 
tion structures of pressure- supported gas clouds in dark 
matter minihalos and lately IGnat fc Sternberg! (|2009l ) con- 
sidered the ionic composition of gas in post-shock cooling 
layers behind fast radiative shock waves. More recently, 
three-dimensional nonequilibriu m modelling of the interstel- 
lar med uim were performed by Ide Avillez fc BreitschwerdtJ 
l|2010al lbT). They considered the evolution of the turbulent 
magnetized interstellar meduim taking into account the 
time-dependent ionization of ten main chemical elements. 
Here, we study how nonequilibrium affects the ionic compo- 
sition and cooling rates of an enriched photoionized gas and 
analyze a parameter space (i.e., temperature, density, metal- 
licity, and ionizing radiation flux) where the nonequilibrium 
effects should be taken into account. 

In the last decade studies of the intergalactic gas be- 
came very intensive. Among the most interesting are the 
conclusions about the tot al mass of baryons i n the Uni- 
verse and its metal fraction. iFukugita fc Peebles! (|2004T l have 
estimated the cosmic baryon budget and found that less 
than 10% of the total baryonic content in the Universe 
is locked in the form of collapsed objects (stars , galaxies, 
groups). More recently iDanforth fc Shulll (120081 s ) founded 
that approximately ~ 50% of the baryons remain unac- 
counted for. A significant discrepancy is also established 



for the total mass of metals observed in the Universe and 
the metal budget estimated from the total star formation 
rate. This la ck of metals is known as the missing met- 
als problem |Pettini 19991 ). Numerical simulations led by 
ICen fc Ostrikei ( 1999h suggest that approximately 30 to 
50% of the cosmic baryons at z = are in the form of a 
diffuse intergalactic medium with temperature of 10 K< 
T < 10 7 K, which has been called a warm-hot intergalactic 
medium (WHIGM ). Subsequent numerical simulations (e.g., 
iDave et al. 2001 ) and more rec ent observational data inter- 
pretation ( Ferrara et al. I [20051 ) also generally support this 
picture. Howev er, a significant part of t he baryons is still un- 
accounted for (|Danforth fc Shulll 120081 ). Thus, the question 
of missing baryons/me tals physical condition is still open 
l|Bregman et al. 1 120091 ') . Many current and future UV/X- 
ray missions (HSlH ChandrC], XMM-Newton__, GALEX0, 
FUSEE World Space Observatory - WSO: SPECTRUM- 
U\fl SPECTRUM-X-GAMMA0, as well as XEUS^) aim 
to find solution of it. In many numerical effor ts for solving 
this question the equilibrium cooling rates (e.g., Cen fc Fand 



l2006t lYoshikawa fc Sasakl200rj ; ICrain et al. 120091 ) are us"ecL 
Obviously, feedback from galactic winds, supernova explo- 
sions and ionizing radiation from galaxies and quasars move 
the WHIGM away from the equilibrium. In this context 
nonequilibrium cooling rates and ionization states are a nec- 
essary component of the WHIGM ionic composition study. 

The paper is organized as follows. In Section 2 we de- 
scribe the details of the model and tests. In Section 3 we 
present our results and discuss possible applications. In Sec- 
tion 4 we summarize our results. 



2 ATOMIC DATA AND MODEL 
DESCRIPTION 

We study the ionization and thermal evolution of a la- 
grangian element of cooling gas. A gas parcel is assumed to 
be optically thin for external ionizing radiation. In our cal- 
culations we consider all ionization states of the elements H, 
He, C, N, O, Ne, Mg, Si and Fe. We take into account the 
following major processes: photoionization, multi-electron 
Auger ionization process, collisional ionization, radiative 
and dielectronic recombination as well as charge transfer 
in collisions with hydrogen and helium atoms and ions. If 
we define the fraction of the ion stage i as Xi, then the 
time-dependent evolution of Xi is determined by 

dXi 
~~dt 



— —fiXi — aiXiUe — /3iXin e 

dt 

f . >He . H . He \ \r 

-{Q nm + Q nHei + % nun + Vi ftHeiijAj 
+a i+ iX i+ in e + fii-\Xi-\n e 

"KCiil^HI + Ci?i n Hcl)A!'i_i 
+ (?7i+l™HII + rjf+iUllel^Xi + l 



(1) 



http:/ /www. stsci.edu/hst/ 
http: / / chandra. harvard.edu/ 
http:/ /xmm. vilspa.esa.es/ 
http: / /www. galex.caltech.edu/ 
http:/ /fuse. pha.jhu.edu/ 
http:/ /wso. inasan.ru/ 

http://hea.iki.rssi.ru/SXG/YAMAL/project_eng.htm 
http:/ /www.esa.int/science/xeus 
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here 7 is the photoionization rate of each ion, given by 



7i = 47T 



hv 



(2) 



and Vifi is the ionization frequency, J„ is the back- 
ground radiation flux, nm, nun, fiHei, WHcii and n e are 
the number densities for neutral hydrogen, ionized hydro- 
gen, neutral helium, ionized helium, and electrons, respec- 
tively , <jj, v is the p hotoi on ization cross sectio n adop ted 
from IVerner et al~l (| 19961 ). IVerner fc Yakovlevl \ 19951 ), a 
is the recombination rate, including radiative recombi- 
natio n (IShull fc van Steenberdll982l; lArnaud fc Rothenflud 



1 19921 ; IVerner et al. I 1 19961; IBadnelll l2006aT). and dielec 



troni c recombination {B adncll 2006b; Bautista fc Badnelll 
20071; IColgan et al 1 120031. 120041 ; IXitun et al. I l20ol. 120061; 



Zatsarinny ct al 
20041; iMazzotta et al 



20031. l2004allbl. 120061; iMitnik fc Badnell 
19981 ; IBadnelll l2006d). and is th e 
collisional ionization rate adopted from Ivoronovl (|l997f ). 

V^j J ? He j C H > C H ° are the charge transfer rates of ioniza- 
tion and recombination with hydr ogen and helium adopted 
from Arnaud fc Rothenflud <|l985l ) and iKingdon fc Ferlandl 



|l996 ), 7j->i is the total photoionization rate including that 



from deep shells followed by Auger-elec tron ejection, the 



Mewd 



(3) 



Auger effect probabilities are taken from iKaastra 
<| 19931 ). 

The system (1) must be complemented by the temper- 
ature equation. Neglecting the change of number of parti- 
cles in the system (for fully ionized hydrogen and helium 
it remains approximately constant) the gas temperature is 
determined by 

dT nY — n e nH^ 
dt Anks 
where n, n e and uh are the electron and total hydrogen 
number densities, T(xi, T, Z, J v ), A(n,T, Z) are heating and 
cooling rates, A is a constant equal to 3/2 for isochoric 
and 5/2 for isobaric cooling, ks is the Boltzman constant. 
Cooling is isobaric when the cooling time is much greater 
then the dynamical time of a system, t c /td S> 1 and it is 
isochoric when the time ratio is opposite, t c /td <C 1. In 
the lack of a dynamic determination we are able to con- 
strain the dynamical time only by the age of the universe, 
t H ~ 2/3H(t). Hot gas (T > 10 6 K) with low overdensities, 
S = (p — (pb(z))) / (pb(z)) ~ 1, has cooling time higher than 
the local Hubble time. In this case the adiabatic approxima- 
tion is a good choice. 

The total cooling and heating rates are calculated 
using the ph o toion ization code CLOUDY (ver. 08.00, 
iFerland et al. I Il998l ). More specifically, we input into 
CLOUDY code a given set of all ionic fractions Xi calculated 
at temperature T, gas number density n and external ion- 
ization flux J(v) and obtain the corresponding cooling and 
heating rates. For t he solar metallic it y we adopt the abun- 
dances reported by lAsplund et al. ] 120051 ), except Ne for 
which the enhanced abundance is adopted (|Drake fc Testal 
l2005h . The solar abundances are listed in Table 1. In all our 
calculations we assume the helium mass fraction Yko = 0.24. 

We are interested in the ionization and thermal evolu- 
tion of gas exposed to external ionizing radiation. A shape 
of ionizing background spectrum gives information about 
population of sources which give major contribution to the 
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Figure 1. The ionizing backgrounds. The Haardt & Madau 
(2001) spectra at redshifts z = 0, 3, 6 are depicted by thick solid, 
dashed and dotted lines correspondingly. The power-law profiles 
with J21 = 1 and the spectral index a = 1.7, 5 are presented by 
thin solid lines. 
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Figure 2. The parameter Un (U is the ionization parameter) 
for the Haardt & Madau (2001) spectra (HM) - solid line, and 
the power-law spectra (PL) - dash lines. 



background. An exposed gas can be either located both near 
galaxy or quasar, or far from any source. In the former the 
ionizing radiation is inherent to a particular source, in the 
latter the gas is radiated by the extragalactic background 
formed by the galactic and quasar population. 

Here we consider two types of spectra: a power-law 
spectrum, which is typical to stellar populations of galax- 
ies, active galactic nuclei and quasars, and the extragalac- 
tic background calculated by usi ng the radiative trans- 
fer in the cosmological simulations l|Haardt fc Madaulll99ej , 
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Table 1. Solar elemental abundances 



Element 


(X/H) 




Carbon 


2.45 


X 


io- 


-4 


Nitrogen 


6.03 


X 


1(T 


-5 


Oxygen 


4.57 


X 


10" 


-4 


Neon 


1.95 


X 


io- 


-4 


Magnesium 


3.39 


X 


io- 


-5 


Silicon 


3.24 


X 


io- 


-5 


Iron 


2.82 


X 


io- 


-5 



l200ll . HM01). A power-law (PL) spectrum shape is J = 
Ji\(y lvH)~ a , where J21 is the flux value at the hydrogen 
ionization limit in units 10 -21 erg s _1 cm~ 2 str -1 , a is the 
spectral index. To study the photoio nization of gas aroun d 
a quasar a is taken equal to 1.7, jSazonov et al. 1 [2004 ). 
but we note that other studies have found both softer and 
harder spectra, with a significant variance about this valu e 
jZheng et al. Ill997l ; iTelfer et al.1l2002t IScott et al. 1|2004 ). 
We note that the quasar spe ctra show a gradual fl attening 
in the 0.2-2 keV band (e.g. ISazonov et al. Il2004) . Such a 
property can be clearly seen for the cumulatuve extragalac- 
tic background for E £ 1 keV at redshifts and 3, where the 
contribution from quasars is dominant. We have also made 
calculations for shallower spect rum, a = 1, because some 
authors use this value (see e.g. iKitavama et al. 1 l200lh . To 
mimic the stellar cluster ionization background, we adopt 
the spectrum with a = 5. For such a spectrum shape, the 
photon-number weighted average intensity above 13.6 eV 
is close to that for the blackbod y spectrum with an e f- 
fective temperature ~ 2 x 10 4 K (|Kitavama et al.l l200ll ). 
Certainly, to model an adequate stellar cluster spectrum 
the theoretical starburst calculations are necessary. How- 
ever, such models depend on many parameters, e.g., initial 
mass function, star formation efficiency, etc. Therefore, we 
have adopted a simple power-law shape. To model the in- 
flue nce from the extr a galact ic ionizing background we use 
the lHaardt fc Madaul (j200lh ionizing radiation flux, which 
is included in CLOUDY code. In our simulations we take the 
radiation spectrum from 1 to 10 4 eV. Figure [1] presents the 
ionizing background profiles used in our calculations. It is 
worth noting that the original spectra presented by HM01 
do not include the absorption in the He II resonant fines 
(|Madau fc HaardtJ l2000l ; iFaucher-Giguere et al~ll2009l ). al- 
though such an absorption can change the abundance s of 
metal ions like CIII and SilV (|Madau fc Haardtl [20091 ') at 
redshifts before the helium reionization. 

To characterize an ionizing background it is convinient 
to use the ionization parameter U defined as the dimen- 
sionless ratio of hydrogen ionizing photons to the total hy- 
drogen density. For a power-law spectrum with the spec- 
tral index a, the ionization parameter can be easy found: 
U = 6.8 x 10~ J J2i /not. Figure [2] shows the parameter Un, 
where n is the total number density, for the extragalactic 
ionizing background (HM01). 

We solve a set of 96 coupled equations (95 for ionization 
states and one for temperature) using a Variable-coefficien t 
Ordinary Differential Equation solver (|Brown et al. Ifl989t ). 
The time step is chosen equal to the minimum value between 




Figure 3. The collisional equilibrium carbon state frac- 
tions for our code (cr osses) are superposed on the results by 
iMazzotta et al.1 lll998T ) (solid lines) and the standard collisional 
equilibrium test of CLOUDY (dash-dotted linos). 
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Figure 4. The collisional time-dependent carbon state fractions 
for our code (thick dashe d lines) are superposed on the results by 
iGnat fe Sternber s (2007) shown by thin solid lines for metallic- 
ities 10 — 3 Zq (upper panel) and solar value (middle panel). The 
roman numbers indicate the ionizaion state. The collisional time- 
dependent cooling rates obta i ned fr om our code (thick dashed 
lines) and IGnat fc St ernberg ( 20071) data (thin solid lines) for 
10 — 3 Zq (lower curves) and solar (upper curves) metallicities are 
shown in the lower panel. 
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(O.lti on , 0.01t coo ;), where t; on and t coo i are the ionization 
and cooling times, correspondingly. 

The evolution of gas certainly depends on a choice of 
the initial condition of gas, namely, the initial ionic composi- 
tion and temperature. We consider su ch a dependence else - 
where (some analysis can be found in lVasiliev et al. ||2010h . 
In this paper we start our calculations from temperature 
T = 10 s K with the collisional equilibrium ion abundances. 
At the beginning we present several tests of the method of 
calculation. In these tests we also start our calculations from 
temperature T — 10 s K with the collisional equilibrium ion 
abundances. During the calculation the gas number den- 
sity is assumed to remain unchanged due to any dynamical 
or chemical processes. We stop our calculations when the 
temperature reaches 10 4 K or photoionization heating rate 
differs from cooling rate by less than 5% (|T — A|/A < 0.05). 
The ionic composition is almost "frozen" at this moment. 
Note that with such a stopping criterion a gas still remains 
"time-dependent" at the end of the calculation, as it has not 
yet reached ionization equilibrium. 

Figures [3]|4] present the main tests of our code used for 
the calculations. At first we tested our model by comparing 
calculations in the collisional equilibrium (we run the simu- 
lations long enough at constant temperature so that eventu- 
ally all ion species settle in ionization equilibrium). Figure [3] 
shows the carbon ionization fractions in the collisional equi- 
librium obtain ed from our cod e (dep icted by crosses) and 
the results by iMazzotta et al. I |l998|) as well as the stan- 
dard CIE test of the CLOUDY code. One can see a good 
coincidence between the results of our code and those ob- 
tained in the previous works. 

Two upper panels of Figure [4] presents the carbon ion- 
ization fractions in the collisional time- dependent case ob- 
tained from our code and the results by I Gnat fe Sternberg! 
(|2007l ). One can see that the fractions in both calcula- 
tions are close to each other for low metallicity value, 
Z — W~ 3 Zq, while the strong deviations can be found for 
solar metallicity. We expect that this difference is due to 
various atomic data and, as a consequence, ionizatio n and 
thermal evolution. Actually, I Gnat fc Sternberg! (|2007r ) have 
utilized the CLOUDY code ver.06.00, whereas in our cal- 
culations more recent version (08.00) of the code is used. 
There were several revisions of atomic data between these 
version^. Also we don't include in our code the cooling rate 
by sulfur, but small variation in the total cooling rate can 
produce different ionic composition. Thus, in low metallic- 
ity gas the ionic compositi o n obt ained in our code is close to 
that of I Gnat fc Sternberg) (|2007h . because the thermal evo- 
lution is determined by hydrogen and helium. Whereas the 
increase of metallicity leads to s tronger deviation between 
our and I Gnat fc Sternberg) ()2007r ) results. Note that the to- 
tal cooling rates in both calculations are close, the rates are 
shown in the lower panel of Figure [4] 



3 RESULTS 

In this section we consider the isochoric cooling rates and the 
ionization states. We leave out of scope of this paper the iso- 

9 http://wiki.nublado.org/wiki/RcvisionHistory 



baric cooling, because the increase of gas density during iso- 
baric process, obviously, leads to the collisionally controlled 
ionization states and radiative cooling (the increase factor 
of density is 10 4 if we study the gas evolution from 10 8 to 
10 4 K). 

3.1 Cooling rates 

In Table 2 we list the main parameters of the models con- 
sidered her J"l. We start from temperature Tj = 10 s K with 
the collisional equilibrium ion abundances. This seems to 
be a good approximation, because the ionic composition at 
so high temperature is mainly determined by frequent colli- 
sions. 

3.1.1 Dependence on metallicity for power-law spectrum 

Figure [S] shows the cooling rates for gas exposed to exter- 
nal radiation for several models listed in Table 2. In the 
panel (a) the dependence of cooling efficiency on the gas 
metallicity is presented. The other parameters are fixed: 
n = 10 -4 cnr 3 , J21 = 1 and a — 1.7. In addition, the 
collisional nonequilibrium cooling rates for the same metal- 
licities are depicted by thin lines. A substantial difference 
between the cooling rates of photoionized and collisional gas 
at temperature lower 10 6 K is clearly seen. The cooling of 
photoionized gas is suppressed in contrast to that of colli- 
sionally ionized gas, because the gas is overionized due to 
photoionization. Indeed, one can see that the cooling rates 
of photoionized gas at T ~ 10 5 K are smaller by an order 
of magnitude in comparison with that of collisionally ion- 
ized gas. An increase of metallicity provides more effective 
cooling of gas and shifts the maximum in the cooling rate to 
higher temperatures. In general, for the same physical state 
of gas the photoionization cooling rates are lower than the 
collisional ones. 

3.1.2 Dependence on flux amplitude and spectral index 

In the panel (b) one can see the cooling rates for gas exposed 
to radiation with fluxes J21 = 0.01, 0.1, 1, the other param- 
eters, a = 1.7, n = 10~ 4 cm~ 3 , Z = Zq, were fixed. The 
cooling rate is getting suppressed significantly with increase 
of the flux. The panel (c) presents the cooling rates for three 
different values of the spectral index a = 5, 1.7 and 1. In 
this case the picture is similar to the previous one. For less 
steep spectrum the cooling rate is suppressed stronger. 

3.1.3 Dependence on redshift for the extragalactic 
spectrum 

The spectrum formed by a cumulative extragalactic radia- 
tion field from galaxies and quasars differs significantly from 
a simple power-law spectrum (see Figure [1}. Consequently, 
the cooling rate of gas exposed to the extragalactic radiation 
also differs from that exposed to a power-law spectrum. Fig- 
ure[5ji presents cooling rates for gas with n — 10~ 4 cm -3 and 

10 The cooling rates and ionic composition for several models 
are available at http://ism.rsu.ru/cool/cool.html 
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Table 2. List of models 



n, cm 3 Z/Zq spectrum 

10" 4 , 1CT 3 , 1CT 2 10~ 3 -1 power-law (PL), a = 1, 1.7, 5; J = (10" 4 -l)J 2 i 
1CT 4 , 10" 3 , 10" 2 10" 3 -1 Haardt & Madau (2001) data (HM), z = - 6 




45674 5 6 74 5 6 7 

log T, K log T, K log T, K 



Figure 5. The cooling rates for gas exposed to ionizing radiation: (a) the dependence on metallicity for gas density n = 10" 4 cm" 3 of 
gas exposed to the power-law (PL) external radiation with J21 = 1 and a = 1.7, the cooling rates for metallicity Z = 10" 3 , 10" 1 and 1 
Zq are shown by thick solid, dashed and dotted lines, correspondingly, the cooling rates for pure collisional case are depicted by thin lines 
of the same type, (b) the dependence on flux amplitude J21 for power-law spectrum with a = 1.7 and gas with n = 10 -4 cm" 3 . Z = Zq, 
the cooling rates for J21 = 0.01, 0.1, 1 are shown by thick dotted, dashed and solid lines, correspondingly; (c) the dependence on spectral 
index a for gas with n = 10~ 4 cm -3 , Z = Zq exposed to J21 = 1, the cooling rate for the spectral index a = 1 is shown by thick solid 
line, the rate for a = 1.7 - thic k dashed line, and the ra te for a = 5 - thick dotted line; (d) the rates for gas with n = 10 -4 cm" 3 and 
solar metallicity exposed to the Haardt fc Madaul l l200ll) background radiation (HM), the solid, dashed and dotted lines correspond to 
the backgrounds at rcdshifts z = 0,3,6; (e) the rates fo r gas with solar metallic ity and number densities n = 10" 4 , 10" 2 cm" 3 (thick 
solid and dashed lines, correspondingly) exposed to the lHaardt &: M adau (200l]) background radiation at z = 3; (f) the nonequilibrium 
(thick) and photocquilibrium (thin lines) cooling r ates for gas with n = 10 ~ 4 cm -3 and metallicities Z = 10" 3 , 10" 1 , 1 Zq (solid, 
dash and dot lines, correspondi ngly) exposed to t he [ Haardt fc M adau (200l)) background radiation at z = 3, the photoequilibrium rates 
are similar to those obtained bv lWiersma et ah ] fe009l) , see details in the text; In all panels, except(a) and (f), the collisionally induced 
cooling rate for gas with n = 10 4 cm" 3 , Z = Zq is depicted by thin dotted line. 



solar metallicity exposed to the HM01 extragalactic back- 
ground spectra at redshifts z — 0, 3, 6. Among these red- 
shifts the ionizing radiation flux is the strongest at z = 3, 
so a difference from the collisional cooling rate is the most 
significant at z = 3. The cooling rates for two other redshifts 
presented here coincide between each other. In comparison 
with the collisional case the rates demonstrate notable dif- 
ferences at T < 10 6 K for z = 3 and T < 8 x 10 4 K for z = 
and 6. Such differences are caused by the overionization of 
the dominant coolants. A more detailed discussion can be 
found in the next subsection. 



3.1.4 Dependence on density and metallicity for the 
extragalactic spectrum at z = 3 

Fugure [5^ shows the dependence of the cooling rates on 
gas density. This dependence is determined by the colli- 
sional ionization and photoionization time ratio, which is 
proportional to n e /J2i- A higher ionizing flux leads to a 
predominance of the photoionization, whereas an increase 
of density intensifies the collisional ionization. For instance, 
the cooling rates of collisional and photoionized gas with 
n — 10 -2 cm -3 (panel e) are equal to each other at 
T Si 5 x 10 4 K. A small difference at lower temperature is 
determined by the overionization of hydrogen, carbon and 




Figure 6. The cooling rates in photo-noncquilibrium (left panels) and photo-equilibrium (right panels) for gas with solar mctallicity 
and n = 10~ 4 cm~ 3 (upper panels) and n = 10~ 2 cm~ 3 (lower panels) , exposed to the extragalactic HM ionizing radiation at rcdshift 
z = 3. The thick solid black line shows the total cooling rate, the thin dash line shows the contribution from hydrogen and helium, the 
other lines (see the legend on the panels) show the contributions from each chemical element. 



oxygen, which are the main contributors to the cooling rate 
in this temperature range. 

Figure [5f shows the dependence on metallicity of 
nonequilibrium cooling rates for gas with number density 
n — 1CP 4 cm~ 3 exposed to the HM01 spectrum at z = 3. 
They differ considerably from the nonequilibrium collisional 
rates presented on the panel (a). The presence of UV radi- 
ation leads to overionized ionic composition, consequently, 
suppresses the cooling rates. 

Recently W lcrsma et al. I l|2009h have calculated the 
cooling rates in the photoionization equilibrium. Because of 
the restricted number of the chemical elements and its dif- 
ferent abundances in our model we re c alcula te photoequilib- 
rium rates similar to lWiersma et al. I (|2009h and compare it 
to the nonequilibrium photoionization cooling rates. We add 
the corresponding cooling rates in Figure [5f- One can see a 
good coincidence between the equilibrium and the nonequi- 
librium cooling rates of photoionized gas for Z ^ Q.IZq. 
This coincidence is a result of the influence of the ioniz- 
ing radiation, whereas for the collisional case the equilib- 
rium and the nonequilibrium cool ing rates differ significantly 
(e.g.. [Sutherland fc Dopita 1993). A sufficiently high ioniz- 



ing radiation flux forces the ionic composition to settle onto 
the equilibrium, though a small difference is seen for solar 
metallicity at T ~ 2 x 10 5 K. Therefore, one can expect 
that nonequilibrium effects will become more significant for 
higher metallicity. 

Let us consider the difference between the nonequilib- 
rium (time-dependent) and photoequilibrium cooling rates. 
Figure [5] shows the individual contributions from each chem- 
ical element to the total cooling rate for gas exposed to 
the HM01 extragalactic ionizing radiation at z = 3 in 
the photo-nonequlibrium (left panels) and photoequilibrium 
(right panels). Figure[6]presents the rates for solar metallicty 
and two number densities, 10 -4 and 10 -2 cm~ 3 . The time- 
dependent rates coincide with the photoequilibrium ones at 
T > 10 6 K. At the same time, for T < 10 6 K the total cool- 
ing rates for both densities differ significantly in the pho- 
toequilibrium and time-dependent models. This difference 
comes from individual contributions of the considered ele- 
ments. One can see that the contributions from C, N, O 
and Ne in the time-dependent case are smaller than those 
in the equilibrium. These elements are overionized in the 
time-dependent case. Such a picture is determined by the 
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Figure 7. The ratio of noncquilibrium to equilibrium cooling rates for gas exposed to a power- law spectrum: a) n = 10 4 cm 3 , 
Z = 1CT 3 Zq and the index a = 1.7; b) n = 1CT 4 cm" 3 , Z = Z G and a = 1.7; c) n = 1CT 2 cm~ 3 , Z = Z and a = 1.7; d) 
n = 10 -4 cm -3 , Z = Zq and a = 5. The ionization parameter U corresponded to the spectrum and number density is shown by extra 
y-axis. The letter "c" on the y-axis (or the area is separated by dash line) corresponds to the ratio of the cooling rates of collisional 
plasma (no photoionization). The shaded area corresponds to the conditions under which the photoheating exceeds the cooling in the 
noncquilibrium case. 



recombination timescales of the C, N, O and Ne ionic states, 
which lead to the recombination lag in the time-dependent 
model. For n = 10 -2 cm -3 the hydrogen and helium should 
be added to the list of the overionized ions. Panels (c-d) 
show a strong difference in the H&He conributions. For 
n = 1CP 2 cm" 3 collisions dominate over photoionization in 
a wide temperature range, so the photo-nonequilibrium cool- 
ing rate coincides with the collisional one at T > 5 x 10 4 K 
(see Figure [5£) and the photoequilibrium and CIE cool- 
ing rates equal to each other at T ^ 4 x 10 4 K. Below 
T~4x 10 4 K the difference between photoequilibrium and 
CIE cooling rates is mainly determined by the overionization 
of hydrogen and it is maximal at T ~ 2 x 10 4 K (the ratio 
of the rates in photoequilibrium to that in the CIE reaches 
a factor of 5). In general, further increase in the number 
density leads to the domimation of collisions in wider tem- 
perature range. It is interesting that for n = 10~ 2 cm -3 
the contributions from carbon and nitrogen are higher than 
these for n = 10~ 4 cm" 3 , e.g. for carbon the increasing fac- 



tor reaches about a factor of 5 around T ~ 7 x 10 5 K, where 
it provides the maximum contribution. 

3.1.5 Importance of nonequilibrium effects 

In order to determine a parameter space, where nonequilib- 
rium effects play a significant role we calculate two sets of 
models: one is for power-law spectra, a = 1.7 and 5, the 
other is for the HM01 extragalactic background. 

Figure [7] presents the ratio of the nonequilibrium to 
equilibrium cooling rates for gas exposed to a power-law 
spectrum with a — 1.7 (panels (a-c)) and a = 5 (panel (d)). 
To compare easily the ratio of the cooling rates of collisional 
and photoionized plasma the ratio of collisional gas (no pho- 
toionization) is shown by the letter "c" on the y-axis (the 
area is separated by dash line) on each panel of Figure [7] In 
the panels (a) and (b) the ratio is shown for Z = 10 _3 Zq and 
the solar metallicities. As one expects the photoequilibrium 
and time-dependent cooling rates are almost equal for the 
lower metallicity. Only for both J 2 i Z 10" 3 and T < 10 5 K 




Figure 8. The cooling rates, carbon and oxygen ionization states in photo-nonequilibrium for solar metallicity and n = 10~ 4 cm -3 
exposed to the power-law ionizing radiation with J21 = 1 (left panels) and J21 = 0.01 (right panels), the spectral index a equals 1.7. 
The thick solid black line shows the total cooling rate, the thin dash line shows the contribution from hydrogen and helium, the other 
lines (see the legend on panels) show the contributions from each chemical clement. 



the rates demonstrate a difference, which originates from the 
efficient H and He recombination in the photoequlibrium 
for such a low ionizing flux. For solar metallicity (panel 
b) the deviation between the rates can be found in wider 
range of both flux and temperature. The time-dependent 
rates tend to the photoequlibrium at J21 <^ f (e.g. the ratio 
is close to 1) that corresponds to the ionization parame- 
ter log U — 0.4 (see the extra y-axis). The increase of 
density, n = fCP 2 cm~ 3 , leads to further expansion of the 
range (panel c), where the photoequilibrium and photo- 
time-dependent rates show a significant difference. Under 
such conditions the photo-nonequilibrium rate tends to that 
in the collisinally cooling gas (the region below the dashed 
line). So the maximum difference between photoequilibrim 
and photo-nonequilibrium rates corresponds to that between 
the CIE and collisinally time-dependent cooling rates (note 
that a significant difference between the collisional equilib- 
rium and nonequlibrium rates at T Si 2 x f 4 K is deter- 
mined by the efficient hydrogen recombination in the CIE). 
Similar difference can be seen in the panel (d), where the ra- 
tio is presented for the steeper spectrum, a = 5. Note that 



in the panels (b) and (d) the ratios are presented for the 
different power-law index a, so the difference between the 
ratios at a given U for the power-law spectra with a = 1.7 
and 5 is due to the different photon distribution (in general, 
the dependence of cooling rates on ionization parameter can 
be considered if the same photon distribution is used). This 
difference is clearly seen for log U = —1.4 and can be ex- 
plained by stronger overionization of high ionic states by 
hard (far UV and X-ray) photons and longer recombination 
lag for the spectrum with a — 1.7 in comparison with that 
for the steepper spectrum with a = 5. The number of such 
hard photons is larger for a = 1.7 than for a — 5, although 
the number of the hydrogen ionizing photons is the same for 
both spectra. 

Panels (a-b) in Figure [8] present the individual contri- 
butions from each chemical element into the total cooling 
rate for gas with solar metallicity exposed to the power- 
law ionizing radiation with J21 = 1 (a) and J21 = 0.01 
(b), the spectral index a equals 1.7. At T > 2 x 10 6 K the 
total cooling rates are very close to each other as well as 
almost equal to the photoequilibrium cooling rates (see Fig- 
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ure[7p). In this temperature range the H, He and Fe ions 
provide the main contribution to the rate. At such a high 
temperature contribution of H and He is due to thermal 
bremsstrahlung. Below this temperature the cooling by the 
Fe, Ne, O and C ions dominate. The contributions from Fe 
and Ne are significant in the ranges T < (0.6 - 2) x 10 6 and 
T (1 — 6) x 10 5 , correspondingly, for both fluxes considered 
here. If the flux (the ionization parameter U, in general) in- 
creases, then the Ne cooling bump strongly diminishes. This 
trend can be clearly seen in the panels (a) and (b). Further 
increase of the ionizing flux leads to the Fe bump diminish- 
ing. In Figure [5] similar situation can be seen for the O and 
espesially for C cooling bumps at T ~ 10 K. Actually, for 
J21 = 1 the contribution from oxygen dominates only at T 
below 10 5 K, but for J21 = 0.01 oxygen becomes a major 
coolant at T ^ 2 x 10 5 . Carbon demonstrates more remark- 
able increase of the cooling rate: for J21 = 1 its contribution 
is negligible, but for J21 = 0.01 it gives a significant part 
of the cooling rate and even becomes one of the dominant 
coolant at T ~ (2 - 4) x 10 4 K. We note that for J 21 = 1 
the difference of cooling rate between photoequilibrium and 
nonequilibrium is low as ~ 10%, so the ionic composition 
is expected to be close to the photoequilibrium: the Ne and 
Fe ions are slightly overionized. The opposite situation is 
for J21 = 0.01, the difference of the cooling rates reaches a 
factor of two at T ~ 2 x 10 5 K. Under such a flux the ions 
of the main coolants (C, O, Ne and Fe) are strongly overi- 
onized in the range T 10 6 K, where the difference of the 
cooling rates becomes significant (see Figure [7p). The Mg 
and Si give a minor contribution in the whole temperature 
range for two fluxes considered. 

Figure HJc-f) present the carbon and oxygen ionization 
states in a gas exposed to the power-law ionizing radiation 
with J21 = 1 (left panels) and J21 = 0.01 (right panels) with 
the spectral index a — 1.7. For J21 = 1 the increase of the 
cooling by oxygen at T ^ 5 x 10 5 K is certainly determined 
by OVI and OV ions (see panel (e)). The contribution from 
oxygen at T 5 x 10 4 K is probably to go down due to 
the existence of the high ionic states only, OV-OVII, whose 
cooling rate is inefficient in this temperature range. Carbon 
is strongly ionized and the CV-CVII ions do not provide 
any significant energy losses. For the lower flux, J21 = 0.01, 
the photoionization rate is smaller and the recombination 
becomes more efficient. The contribution from carbon at 
T £ 2 x 10 4 K is determined by the rapid increase of the 
CIII and CIV ion abundances. The dominant contribution 
from oxygen in the wide temperature range can be explained 
by the OV-OVI ions at T > 10 s K and the OIII-OIV ions 
at T ^ 10 s K. Note that for a photoionized gas the decrease 
of metallicity leads to increase of the H and He contribu- 
tion to the cooling rate. The cooling of low metallicity gas, 
Z ~ lO _3 i?0, controlled by H and He ions only, the contri- 
bution from other elements is negligible. 

Figure [5] demonstrates the ratio for the HM01 back- 
ground spectrum at redshift z. For the extragalactic HM01 
background the picture is similar to that for a power-law 
spectrum: an increase of the metallicity leads to an increase 
of the magnitude of deviations from the equilibrium. For 
low metallicity, Z ~ 10~ 3 Zq, the ratio remains close to 
unity in the wide range of temperature. Only for T ~ 10° K 
it deviates from unity about 20-30% that is determined by 
the overionization of helium (see Figure [5f). Figure [9] shows 
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Figure 9. The ratio of nonequilibrium to equilibrium cooling 
rates for gas with density n = 10~ 4 cm -3 and solar metallicity 
exposed to the Haardt & Madau (2001) background spectrum at 
redshift 2. The shaded area corresponds to the conditions under 
which the photohcating exceeds the cooling in the nonequilibrium 
case. 



the ratio for solar metallicity. As one expects the nonequi- 
librium effects are important for T ^ 10 6 K in the whole 
redshift range, z = — 6. The maximum deviation reaches a 
factor of ~ 2 at T ~ 10 5 K for z ^ 1 and z ^ 5. As the extra- 
galactic background is maximal around redshifts z = 2 — 3 
a saddle point on the surface of the cooling rate ratio is ob- 
served: the deviation from the equlibrium is minimum in this 
redshift range. As we mentioned above the difference from 
photoequilibrium is determined by smaller contribution to 
the cooling from C, N, O and Ne ions, because they are ove- 
rionized. The ionizing flux in z = 2 — 3 is the highest, so 
the ionic composition is moved to photoequlibrium (see the 
previous section). For both 2 ~ and z ~ 6 the ionizing flux 
decreases and the difference from photoequilibrium becomes 
more clear. 

Let us summarize where the photoequilibrium calcula- 
tion can be used and how the parameters (metallicity, tem- 
perature, density and UV intensity) should be changed so 
that the nonequilibrium effects become important. The cal- 
culations presented above show that the photoequilibrium 
is reached in high-temperature gas, T £ 3 x 10 6 K, for any 
metallicity or ionization parameter. Note that the photoe- 
quilibrium equals to the CIE for such a high temperature. 
Below T ~ 3 x 10 6 K the deviation from the photoequi- 
librium grows. For a metallicity as low as 1Q~ 3 Z@ such a 
deviation is minimum. The increase of metallicity leads to 
the widening of the parameter space, where nonequilibrium 
effects are significant (Figures [7] and |5J) , that is clearly seen 
for solar metallicity. This can be easily explained by the 
overionization of metal ions and the increase of the con- 
tribution from metals to the cooling rate for higher metal- 
licity. For collisionally ionized gas the deviation from the 
equlibrium at T < 10 6 K is significant for any metallicity. 
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Table 3. The transiton from photo-EQ to NEQ at T < 10 6 K. 



the conditions, where 
\EQ ^ f^NEQ ; g expectec l 


a variation of parameters 
leads to A E Q ^ jyNEQ 


high ionization parameter U 
(high UV flux/low density), 
U > 10~ 2 for Z < lCr 2 Z , 
U £ 1 for Z ~ Zq 


either increase of Z 

or decrease of U parameter 

(decrease of flux/increase 

of density) 



In a gas exposed to ionizing radiation the photoequilibrium 
is reached faster than higher ionization parameter. At low 
ionization parameter (low flux intensity /high density), the 
differences between equilibrium and nonequilibrium indeed 
become large, but these are just the differences in the colli- 
sional case. Thus, in Table |3] for T < 10 6 K we provide the 
parameter range under which the cooling rates in the equi- 
librium (EQ) and nonequilibrium (NEQ) are close, and a 
change of parameters, which violates this equality. Here we 
should also pay attention on that the ionization parameter 
U does not fully describe the spectrum and in general the 
transition from photoequilibrium to nonequlibrium depends 
on the ionizing spectrum shape. 

As a conclusion we note that the deviation of the cool- 
ing rates between photoequilibrium and nonequilibrium is 
determined by recombiantion lag in the latter case. The 
maximum difference is for collisionally controlled gas. Such 
condition in gas is reached for ionization parameter as low 
as log U — 5 (see Figures [7] [9} ■ The difference decreases 
with the ionization parameter growth. For sufflcietly high 
ionization parameter the photoequlibrium can be reached. 
However, the transition region from photoequilibrium to 
nonequilibrium is quite wide. For ionization parameter as 
low as log U ~ — 5 the nonequilibrium collisional cooling 
rates and ionic composition are a better choice. The pho- 
toequilibrium is reached for ionization parameter as high as 
log U > - 2 for Z < 1Q- 2 Z Q and log U>0 for Z ~ Z Q 
(see Table [3]). For supersolar metallicity the photoequilib- 
rium is expected to reach for higher ionization parameter. 
So where the physical conditions are favour to collisional 
ionization (the collisions dominate or its contrubution to 
the total ionization rate is significant), the nonequilibrium 
(photo) ionization calculations should be conducted. 

3.2 Ionization states of metals and their ratios 

3.2.1 Collisions or photoionization: CIII and CIV ions 

In the previous sections we have considered the time- 
dependent cooling functions of a photoionized gas, their de- 
viations from the photoequilibrium rates and the contribu- 
tion to the cooling rate from chemical elements. We have 
found that the transition region from photoequilibrium to 
nonequilibrium is quite wide. Obviously, the reason of such 
a transition is the recombination lags in the nonequilirium 
model. Consequently, a significant difference of the ionic 
ratios is expected between the photoequilibrium and the 
nonequilibrium. Moreover, which of two processes, collisions 
or photoionization, is dominant can be also identified using 
ionic ratios. For gas with a given density and metallicity 
exposed to an ionizing radiation the transition between col- 
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Figure 10. The ratio of the CIV ion fraction in the nonequilib- 
rium case to that in the equilibrium (left panels), and the ratio, 
R = x(CIII)/a;(CIV), in the nonequilibrium to that in the equi- 
librium (right panels) for gas with density n = 10 -4 cm -3 and 
solar metallicity expo sed to the power-law spe ctrum with a = 1.7 
(upper row) and the lHaardt fc Madaru d200lj) background spec- 
trum at redshift z (lower row). In the panels (a) and (c) The letter 
"c" on the y-axis (or the area is separated by dash line) corre- 
sponds to the ratio of the cooling rates of collisional plasma (no 
photoionization). The shaded area corresponds to the conditions 
under which the the ratio does not exist (photohcating exceeds 
the cooling in the nonequilibrium case) or the ratio is out of the 
ranges presented in the gray scale. 



lisionally controlled or photoionization regimes depends on 
the photoionization rate, which is determined by spectrum 
shape and flux magnitude. 

Figure [TU] shows the ratio of the CIV ion fraction in 
nonequilibrium to that in photoequilibrium (left panels) for 
gas with density n = 10 -4 cm -3 , solar metallicity exposed 
to the power-law spectrum with a = 1.7 (upper) and the 
HM01 background spectrum at redshift z (lower). As one ex- 
pects the CIV fraction in the photoequilibrium strongly dif- 
fers from that in the time-dependent case: the ratio reaches 
several times. The significant deviation from photoequilib- 
rium can be also found for the HM01 background. In the 
panel (a) one can see that the decrease of the UV flux am- 
plitude leads to that the CIV fraction ratio tends to that 
in collisionally controlled gas. We should note once again 
that due to the recombination lag the carbon ionic states 
in the nonequilibrium are always overionized in comparison 
with these in the equilibrium (see the panels (a) and (c)): 
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Figure 11. The dependence Nqiv/Novi versus Njqy /Novi 
in the photo-nonequilibrium model for cooling gas with n = 
10 — 4 cm -3 and Z = Zq exposed to the ionizing radiation with 
a = 1.7 and two fluxes J21 = 0.01 and 0.1 (the corresponding 
labels are shown). The same dependence in the photocquilibrium 
are by triangles connected by solid line for J21 = 0.1 and by 
dashed line for J21 = 0.01. Gas temperature is indicated by gray 
scale along the trajectories: from hot (black) to cold (light gray) 
gas. The data point s show the ionic ratios observed in metal ab- 
sorbers (see Table 8 lFox et al."1l2006l) . 



where the ratio is below 1, the CV state dominates in the 
nonequilibrium, e.g. the CV is overionized, and in the oppo- 
site case, where the ratio is greater than 1, the CIV is ove- 
rionized, whereas in the photoequilibrium under the same 
conditions the CIII fraction increases fastly. This is clearly 
seen in the right panels of Figure 1101 where the ratio of the 
values, R = x(CIII)/a;(CIV), in the nonequilibrium to that 
in the equilibrium is presented. Indeed, in the panel (b) the 
shaded region in the bottom left Conner corresponds to the 
ratio less than 0.1, i.e. the denominator of the ratio, the R 
value in the equlibrium, is high. 

In the panel (d) the strong deviations of the R value 
from the photoequilibrium can be found for the HM01 
spectrum. Because the CIII/CIV ratio is used as one 
of the important observable indicator to discriminate be- 
tween the Hell reionizaton models dSchave et al. I 120031 ; 
lAgafonova et al. 1120071 ; iMadau fc Haardtll2009l ). we should 
pay attention on the significant difference at redshifts 2 = 
2 — 3, when the Hell reionization is expected to occur. The 
nonequilibrium analysis of t he CIII/CIV ratio in the ab- 
sorption spectra obtained by ( Agafonova et al.ll2007l ) can be 
found in (|Vasiliev et al. Il2010h 7 Certainly, for lower metallic- 
ity smaller deviations are expected. However, one can con- 
clude that the nonequilibrium effects are significant and it 
should be taken into account in the simulations of the IGM 
state. 



3.2.2 An example of the ratio of column densities 

The ratio of colu mn densities of ion i of e lement m and ion 
j of element n is (|Gnat fc Sternberdl2007l ) 



where Am, A n are the abundances of elements m and n (rela- 
tive to H), Xi,Xj are the fractions of corresponding ionization 
states. The absence of the hydrogen density dependence pro- 
vides more general tool fo r diagnostics of the IGM. Following 
iGnat fc Sternberg] (|2007l ) we consider the ratios Nciv /Nqvi 
and A r Nv/A r ovi- 

Figure [TT] shows these ratios for power-law spectrum 
with a = 1.7 and two fluxes J21 = 0.1 and 0.01 in photoequi- 
librium and photo-nonequilibrium. The "cooling tracks" are 
the gray diagrams, where the gray gradation corresponds to 
gas temperature. As can be expected, the strong difference 
between the photoequilibrium and time-dependent tracks is 
seen for J21 = 0.01, whereas the tracks for J21 = 0.1 become 
closer. Indeed, the ionic ratios in a gas exposed to strong ion- 
izing radiation significantly differ from the ratios in collision- 
ally controlled gas (see e.g. Figure lTOjl . The decrease of ioniz- 
ing flux increases the role of collisions and leads to the transi- 
tion between photoionized and collisional plasma. However, 
as one can see in Figures 171 and I5IHT1 there is no sharp tran- 
sition between photoionized and collisional plasma. So ionic 
ratios observed in an object can be satisfied by both pho- 
toionized and collisional time-dependent conditions, e.g. for 
J21 = 0.001 the ratio, R = a;(CIII)/a:(CIV), demonstrates a 
negligible difference from the collisional time-dependent ra- 
tio (see Figure ITTJb ). But the difference becomes very clear 
for J21 = 0.01 at T < 10 5 K. 

Figure [TT1 shows the observational data for t he high ve- 
locity clouds in the Galac ti c halo (see Table 8 in lFox et al. I 
120061 ). IGnat fc Sternberg! l|2007l ) showed that the column 
density ra t ios in the high velocity clouds observed by 
iFox et al. I (120061 ) can be fitted using the nonequilibrium 
collisional model. Here we confirm that the nonequilibrium 
photoionization models can also explain the observed ratios, 
for a gas with n ~ 10~ 4 cm~ 3 exposed to J21 = 0.01 — 0.1 
and a = 1.7 in particular. For example, the ionization states 
of carbon and oxygen for J21 = 0.01 can be found in Fig- 
ure [5] Under such a flux the CIV, NV and OVI fractions are 
close to these in the collisional time-dependent model, but 
strongly overionized in comparison with the photoequilium 
ones. More detailed study is out of scope of this paper. 



4 CONCLUSIONS 

In this paper we have presented the nonequilibrium 
cooling rates and ionization states of gas enriched with heavy 
elements and photoionized by external ultraviolet /X-ray ra- 
diation. We have studied the dependence of cooling rates on 
the gas density and metallicity, and also by the radiation 
field (power- law and the extragalac tic background spectra 
obtained by lHaardt fc Madaul 200 ll ) . The calculations are 
done for a wide range of physical parameters: 10 4 — 10 8 K 
in the gas temperature, Z = 10~ 3 — 1 Zq in the metallicity. 
For the power-law spectrum, the spectral index a is taken 
equal to 1, 1.7 and 5, and for the extragalactic background 
spectra the redshift varies from to 6. In all calculations, 
we assume that the gas is optically thin. We have found the 
following. 



The time-dependent cooling rates depend significantly 
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on the gas metallicity and density as well as on the flux and 
shape of the ionizing radiation. 

• The time-dependent cooling rates and ionic composi- 
tion of gas differ strongly from those in the photoequilib- 
rium due to the overionization of ionic states in the time- 
dependent case. The time-dependent cooling efficiency is 
generally lower than that in the photoequilibrium, and the 
nonequilibrium effects can be clearly seen at T < 10 6 K for 
the solar metallicity. 

• The difference between the time-dependent and equi- 
librium rates can be significant, its magnitude is maximal 
at ionization parameter as low as log(7 £ — 5, and similar 
to the difference between the equlibrium and the nonequi- 
librium cooling rates in the collisionally controlled gas (see 
Figures Oil- 

• Contributions by individual chemical elements to the 
total cooling rate can be significantly different for time- 
dependent and equilibrium plasma at T 1 6 K. 

• The mismatch of both ionic states and their ratios in 
the photoequilibrium and photo-nonequilibrium can reach a 
factor of several. 

In this paper we have considered the external radiation and 
gas parameters, under which the cooling rates and ionic com- 
position in the photoequilibrium and photo-nonequilibrium 
are expected to be close to each other, and a change of 
the parameters which violates this equality. The parame- 
ter space where the nonequilibrium effects are significant 
is quite wide. Because the difference between the photoe- 
qulibrium and the nonequilibrium is maximal for the ion- 
ization parameter log U ,$ — 5, then the nonequilibrium 
collisional cooling rates and ionization states are a better 
choice for such a low ionization parameter. We have found 
that the difference decreases with the ionization parame- 
ter growth, and the photoequilibrium is reached for ioniza- 
tion parameter as high as log U £ — 2 for Z £ 10 -2 Zq and 
log U %j for Z ~ Zq (see Table 0J. For supersolar metal- 
licity the photoequilibrium is expected to reach for higher 
ionization parameter. We have concluded that the nonequi- 
librium photoionization calculations should be conducted, 
when the contribution from collisions to the total ioniza- 
tion rate is significant. The main reason of the deviation 
in the time-dependent case is the recombination lag, so the 
ionic states are expected to be overionized, especially this 
is notable for metals. Therefore we should pay attention on 
the importance of use nonequilibrium cooling rates for near- 
solar (and above) metallicity of gas exposed to arbitrary 
ionizing radiation flux. Thus, the nonequilibrium cooling 
rates are the necessary component for study metals in the 
1GM, mixing metal processes, damped Lyman-a systems, 
galactic superwinds and so on. Also to simulate correctly 
the WH1GM ionic composition we need nonequilibrium pro- 
cesses and more exact models for ionizing background flux. 
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